Measuring functional renormalization group fixed-point functions for pinned manifolds 
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Exact numerical minimization of interface energies is used to test the functional renormalization group (FRG) 
analysis for interfaces pinned by quenched disorder The fixed-point function R(u) (the correlator of the coarse- 
grained disorder) is computed. In dimensions D = d + 1, a linear cusp in R"{u) is confirmed for random bond 
(d = 1, 2, 3), random field (d — 0, 2, 3), and periodic (d = 2, 3) disorders. The functional shocks that lead to 
this cusp are seen. Small, but significant, deviations from 1-loop FRG results are compared to 2-loop corrections. 
The cross-correlation for two copies of disorder is compared with a recent FRG study of chaos. 
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Systems with quenched (frozen-in) disorder often exhibit 
glassy phases at low temperature. Standard perturbative meth- 
ods fail to describe these phases and exact results are limited 
to ID and mean field models iii. It has been quite a 
challenge to develop field theoretic and renormalization group 
(RG) methods, which must include both multiple metastable 
states and spatial fluctuations in finite dimensions, to describe 
universal properties of these phases. Proposed field theo- 
ries are unconventional and harder to control than those de- 
veloped for pure critical systems. An expansion around the 
mean-field replica symmetry (and ergodicity) broken (RSB) 
solution, much studied in spin glasses, is very difficult even 
at tiie 1-Ioop level 0]. The functional RG (FRG) was de- 
veloped for elastic objects pinned by substrate disorder and 
random fields. This class has numerous physical realizations, 
including vortex lattices, magnetic systems, and charge den- 
sity waves Jsl 0,0,1^. The 1-Ioop FRG has been extended 
to describ e, e.g ., depinning of a driven interface [9J, activated 
dynamics llOll . quantum models lITlll . and sensitivity of con- 
figurations to disorder changes ("chaos") fl^. Since the FRG 
parameterizes the effective action by functions, rather than the 
few couplings of standard RG, it is better suited to handle an 
infinite number of marginal parameters at the upper critical 
dimension (or runaway flows as in correlated fermions 1 13|). 

When applying the FRG to pinned elastic manifolds param- 
eterized by a scalar displacement field u{x), the function in 
the effective action whose flow is relevant below d = 4 is de- 
noted by R{u). Physically, this function represents a coarse 
graining of the coiTelator of the pinning potential; it encodes 
an infinite number of couplings, _R^^"'(0), n = . . . oo. An 
unusual feature of the theory is that R"{u) can develop a lin- 
ear cusp around m = at finite scale |6]. In the space of 
non-analytic functions, perturbative control was recovered to 
one-loop order (i.e., to 0(e = 4— d)) and fixed-point functions 
R{u) obtained for various universality classes |6, 8, 14]. The 
relations between this cusp singularity, multiple metastable 
states and shocks in energy landscapes have been vividly de- 
scribed 1 15|. The FRG agrees with phenomenological mod- 
els and successfully predicts the roughness exponent ( of the 
pinned interface, with the d isorder-aver aged correlation func- 
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lll4ll whether the FRG can be extended in a systematic loop ex- 
pansion, i.e., to higher order in e. Dealing with a non-analytic 
action is very subtle 1191 . and even 1-Ioop consistency is not 
obvious Recently, a candidate renormalizable field the- 
ory for statics fl^.'iT] (and a distinct one for depinning fl^) 
was obtained beyond one loop. Crucial to its construction is 
the property that the cusp remains linear to higher orders. If 
confirmed, the FRG provides a simpler method to attack glass 
problems where the RSB phenomenology can be avoided. 

This paper presents a new level of "smoking gun" tests of 
the FRG for manifolds, by directly measuring the fixed-point 
function R[u) for three universality classes (Figs. [Tl3t . This 
is achieved, as proposed in Ref. 12311 . by adding to the dis- 
order a parabolic potential (i.e., a mass m) with a variable 
minimum location v. The resulting sample-dependent free 
energy V{v) defines a renormalized potential whose second 
cumulant coiTelator in v space is the same R{v) function as 
defined in the field theory (from the replicated effective ac- 
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tion {u{x) - u(0))2 ~ x^'- ll^llTllli 

Though much has been achieved, it has been questioned 
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FIG. 1: Filled symbols show numerical results for Y{z), a normal- 
ized form of the interface displacement correlator —R"(u) [Eq. 
for D = 2 + 1 random field (RF) and _D = 3 + 1 random bond 
(RB) disorders. These suggest a linear cusp. The inset plots the nu- 
merical derivative Y'{z), with intercept ^'(0) ~ —0.807 from a 
quadratic fit (dashed line). Open symbols plot the cross-correlator 
ratio Ys{z) — Ai2(2:)/Aii(0) between two related copies of RF 
disorder It does not exhibit a cusp. The points are for confining 
wells with width given by AI^ = 0.02. Comparisons to 1-loop FRG 
predictions (curves) are made with no adjustable parameters. 
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tion - deviations arise only in higher cumulants (t^). This is 
analogous to measuring the coupling constant and the distri- 
bution of total magnetization in pure systems, which underlie 
phenomenological RG and finite size scaling |24]. As in pure 
systems, the FRG predictions are universal at coarse grained 
scales, but require specifying the large scale BCs. The mass 
provides these conditions and also allows one to control and 
quantify the zero mode (center of mass) fluctuations, yielding 
the coupling function R{u). The same procedure allowed an 
exact calculation |23| of R{v) for the D = + 1 theory with 
RF disorder (Sinai's model). 

We numerically compute the FRG zero temperature fixed- 
point functions using exact ground state configurations. We 
study interfaces embedded in dimensions D — d + 1, d — 
0,1,2,3, including random bond (RB), random field (RF), 
and periodic (RP) disorder universality classes. We focus 
on universal, parameter free functions; treatment of universal 
amplitudes requires more details and is presented separately 
llisll . The linear cusp in A{u) = —R"{u) is confirmed in all 
cases. For periodic disorder, A(m) is consistent with the con- 
jectured parabolic form. For RB and RF disorder, the scaled 
A(u) are distinct from the 1-loop calculations and are closer 
to the two-loop results, though the curves exhibit at most a 
weak dependence on d. The functional shocks responsible for 
the cusp in A(u) are directly seen. The higher statistics of 
these shocks are consistent with d — Q Burgers intermittency. 
Cross-correlation (chaos) fixed points for two related copies 
of the disorder show a rounding of A(w) that is consistent 
with recent FRG predictions 1 12]. 

The continuum Hamiltonian for an interface u{x) of inter- 
nal area 57 with elastic constant K, confined in a parabolic 
well centered at v is 

(1) 

where the random potential V has correlations 
V{0,x)V(u,x') = Ro{u)5('^\x - x'). The RB univer- 
sality class describes short ranged Ro{u), the RF class has 
Rq{u) ^ —flul at large u, while the RP class describes 
periodic correlations, e.g., Ro{u + 1) = Rq{u). The bare 
correlator Ro{u) becomes R{u) upon coarse graining. Given 
a UV cutoff scale b, for fixed ilb^'^, and continuous V{x, u), 
the minimum energy configuration u{x; v) is unique and 
smoothly varying with v, except for a discrete set of shock 
positions where u{x; v) jumps between degenerate minima. 

For numerics, interfaces u{x) are described by a set / of 
edge-sharing plaquettes p. Plaquettes are dual to the edges in 
a regular lattice composed of H layers. Each layer has 
points, unit cell volume fio, and periodic BCs. Each point is 
connected to points in the layer above by k bonds, so that an 
interface / has kL"^ plaquettes |25]. The energy Tiiatt of /, 
confined by a well of strength Af centered at v, is 

Wlatt(«) = E {^["(P) - + ^(P)| ' (2) 
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where u{p) is the layer index for plaquette p and U{p) is 
the disorder potential. Long-wavelength elasticity arises from 
combinatorial effects 1 16]. For RB disorder, U (p) is a Gaus- 
sian variable h{p) with zero mean and variance dg = 1, while 
for RF disorder U{p) is the sum of h{p) along a path of edges 
connecting p to the bottom layer RP disorder with period 
P is constructed by stacking H/P identical RB samples of 
thickness P. Given U (p), v, and M, the ground state /gg is 
found using a program that accommodates all lattices, dimen- 
sionalities, and disorder types. The new version of the core 
max-flow algorithm 116.] in our code has been directly tested 
against standard Ubraries 12611 and earlier energy minimiza- 
tion calculations |16, 11811 . The height H is always large 
enough that the finite size effects are controlled only by L 
and M. Lattice discreteness is evident at high values of M, 
so we choose M < 0.2. Continuum and discrete models are 
then related by equating energies Tiiatt and H, displacements 
u{p) and u{x), interface areas H.QKb'^L'^ = H., well strengths 
TO^ = M2(f]g5'i)-i, and disorder strengths cr = t^^. The 
effective elastic constant K was also measured |25]. 
We computed the discrete force-force correlation i23ll : 

Auttiv) = M^{kL'')[v' + v~ uoiv' + v)][v' - uoiv% (3) 

where the mean position uq{v) = {kL'^)~^ J2pei "(p)- The 
averaging (overline) is for N > 10^ samples with < v' < P 
for RP disorder; RF and RP samples are self-averaging over 
v' (we slide v' over more than 10^ times the interface width 
while computing minima in a window of thickness _ff w 20 
centered at v'). The plots we present have la error bars com- 
puted using direct resampling of the data and are thus ex- 
pected to overlap the large sample number limit with a prob- 
ability of 68%. To check our procedure, we confirm that 
du Aiatt {u) is consistent within errors with the value o-q 
for RF disorder and with the value for RB and RP disorders. 

The FRG predicts that, for large volumes il/b'^, the rescaled 
correlator A(z), defined by A(u) — m'^^'*^ A(Mm^), con- 
verges as m — > to the FRG fixed-point function A*(z), 
which depends only on d and disorder. Using Eq. Q, conver- 
gence of Af^'^-'^AiattC^M"^) was evident for L > 16 (L > 8 
ind = 3) andM < 0.2, choosing 1 1^] C = 2/3, 0.44, 0.22 for 
d = 1, 2, 3 RB disorder and ( = (A - d)/3 for RF disorder 
The interface widths grow slowly (^ = 0) for RP disorder 
As the FP functions still contain an amplitude and a scale, we 
introduce the normalized function Y(z), 

A{u) - A(o)rKO, (4) 

so that Y{0) = 1 and with scale ^ chosen according to 
disorder type: for the RP model, ^ = P, for RF disor- 
der, ^ is set so that dzY{z) = 1, and for RB disorder, 
Jp°° dzY'^{z) = 1. This function is predicted to be univer- 
sal with a dependence in d, Y{ z\ d), that can be computed to 
second order in e 4 — c? I19ll2lll 

Y{z-d)^Yi{z) + eY2{z) + 0{e^), (5) 
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with Yi{z) the 1-loop estimate f^.'s','!?]. Computation of K 
is required to fix universal information not retained in Y{z), 
e.g., the ampHtude A(0) for RF disorder i25ll . 

We plot illustrative examples of Y{z) in Fig.Q] In all cases, 
an apparently linear cusp is found for Y(z) (with finite inter- 
cepts for fits to F'(z)). The normalized functions are remark- 
ably close to 1-loop predictions, with no adjustable parameter 
We now turn to a systematic analysis of these functions, their 
deviation from 1-loop results, and related data. 

We start with RF disorder. The ERG predictions for the 
functions Yi(z) and 1^2(2) in © obtained from lineariz- 

I relation jEH] 
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^Y-l-\nY-^F{y) 



where 



Jo' dy^y-l-lny-§Fiv) 

F{y) = 2y - 1 + ^ - ^Iny + Li2(l - y). Plots of 
the differences Y{z) — Yj (z) between the numerical result 
and the 1-loop prediction 16], for several sizes and masses in 
D — 2+1 and D = 3+1, are shown in Fig.|2] There are small, 
but statistically significant, systematic deviations from Yi{z). 
The sign of the expected corrections linear in e, Y2 (z), agrees 
with numerics. This function changes sign at Zc — 1.668 . . ., 
near the observed location. The magnitude of eY2{z), setting 
e = 1, is nearly consistent with numerics for all d. We in- 
clude 0+1 numerical results (compatible with Refs. I23IE7I1 ') 
for comparison. Points for D — 2+1,3 + 1 are both close 
to D = + 1 results. Our computed slopes at the origin, 
-Y'iO) = 0.815(7) (3H-1) and -Y'{0) = 0.811(6) (2+1), are 
to be compared with the FRG value 0.7753 . . .+ (0.0328 . . .)e 
and the d = 12311 value 0.8109 . . .. The near equality of the 
d = curve and I2 (z) appears to be a coincidence. Although 
more work is needed to resolve the differences (e.g., d — 
from d = 2, 3) the trend of the FRG results is encouraging. 

For RB disorder, R{u) is expected to decay (so A(u) has a 
zero). Fixing ^ as stated sets a non-universal scale. The differ- 
ences Y{z) — Yi{z) are plotted in Fig.|2| we again find small 
but significant deviations from the 1-loop prediction, with at 
most a weak dependence on d (within error bars). The O(e^) 
expansion in this case is found from series and numerical so- 
lutions 1 21]. The resulting 12(2) again agrees well in sign 
and shape with the data, with a magnitude given by e sa 1 . We 
have constructed 2-loop interpolations which agree with the 
data in all d |25]. The situation resembles that for RF disor- 
der, even though deviations have the opposite sign. 

Results for the function Y{z) for RP disorder are shown in 
Fig.HJfor d ~ 3; similar results hold for d — 2. The 1- and 2- 
loopFRGspredict|8,21] aparabolicform, A(u) = A(0)(1- 
6u{l — u)), as do the d = and the large-d cases (with a single 
shock as TO ^ fisl] and many small independent shocks per 
period |20], respectively). Counting of derivatives in the FRG 
equation has also indicated that the parabolic form holds for 
any finite d. The parabolic form is consistent with our results 
as TO ^ 0. 

The use of a harmonic well allows one to define and study 
the shocks in the force landscape. As v increases, sections of 
the manifold have degenerate minima at positions Vg and the 
polarization jumps forward by J d'^x [u{x;v'^) — u{x]v~)]. 
These are shocks in a functional (scalar for d = 0) de- 
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FIG. 2: The difference between the normalized correlator Y{z) and 
the 1-loop prediction Y\{z) for RF disorder in D = (1, 2) + 1 and 
RB disorder inD = (l,2, 3) + l. The dashed lines are the linear 
2-loop correction Y2{z) — ''"^Jf^ \e=o of Eq. |3. For each disorder 
class, the data are close to each other and to the d = and e = 1 
linear two-loop estimates, but are distinct from the 1-loop result. 



caying Burgers equation |23|], with the renormalized force 
V — u{v) corresponding to velocity and m^^ corresponding 
to time in Burgers turbulence. Examples of these disconti- 
nuities in the renormalized force are shown in the inset of 
Fig. 13 We have seen shocks merge as m decreases j|25j]. 
The moments of the renormalized force are Sn{v — v') = 
{v — v' — u{v) — u{v'))"-. A linear cusp in S2 is confirmed 
by our study of A(u). A prediction of the FRG in d > Q 
is that S^{v — v') ^ {v — v') at small v — v' , in accord with 
exact results for d — Q. Linearity of all Sn, n > 2 is a hall- 
mark of intermittency in d = Burgers turbulence. Our data 
show linearity of 6*3 (Kolmogorov scaling) and S4 in v — v' 
for cases studied. This indicates that shocks do not cluster 
beyond simple statistical fluctuations. 

When the pinning potential is perturbed, correlations be- 
tween the original and perturbed samples remain for RF dis- 
order and are described by a new chaos FRG fixed point 
lll2ll . We test this prediction using related disorders Ui(ji) 
and U2{p) = [Ui{p) + S ■ W{p)]/\/l + S'^, where the per- 
turbation W{p) is a mean-zero univariate Gaussian and 5 is 
the perturbation strength. We measured the cross correla- 
tor Ai2(f; - v') = kL'^M'^{v - uo,i{v)) {v' - u q,2{v'))- We 
check the sum rule Ai2(u) = cr/\/l + and normal - 

izeviai;(z) = Ai2(^z)/Aii(0). We find (Fig.Q that K,(z) 
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FIG. 3: (a) Plots of the normalized pinning force correlator Y{z) for 
RP disorder in d = 3. For these values of m, the period P = 4 
points have mostly converged to the parabolic RP fixed-point func- 
tion, while the P = 8 curves are still crossing over from the RB 
to the RP universality class, (b) Residuals relative to the parabolic 
shape vanish, within error bars, for larger sizes and P = 4. 
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FIG. 4: A plot of the normalized cross-correlator , computed 

for RF disorder in d — 2, 3, showing the sensitivity to disorder of 
magnitude S, compared to the l-loop prediction and to numerical 
D = + 1 computations (error bars not shown; la errors are about 
1/2 of symbol size). Inset: the renormalized pinning forces v—uo (v) 
for a sample (solid line) and aS = 0.1 perturbed sample (dashed line) 
in a typical sample; cross-correlations of such data give the main plot. 



is rounded, as predicted 1 12|. The computed Ys{0) is near the 
l-loop prediction (see Fig.0}. 

Our numerical calculations confirm the main features of the 
FRG approach to the glassy system of pinned interfaces, es- 
pecially the shape of R"{u) and its linear cusp, for a variety 
of disorders and dimensionalities. FRG computations to 2- 
loop order significantly improve upon the l-loop results. The 



functional shocks found are consistent with expectations; their 
statistics merit further study. 
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